Maximum Angle of Stability of a Wet Granular Pile 
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Anyone who has built a sandcastle recognizes 
that the addition of liquid to granular materials 
increases their stability. However, measurements 
of this increased stability often conflict with the- 
ory and with each other 0, 0, IE IE IE IE 01- 
A friction-based Mohr-Coulomb model has been 
developed |E Igj- However, it distinguishes be- 
tween granular friction and inter-particle friction, 
and uses the former without providing a phys- 
ical mechanism. Albert, et al. [2\ analyzed the 
geometric stability of grains on a pile's surface. 
The frictionless model for dry particles is in ex- 
cellent agreement with experiment. But, their 
model for wet grains overestimates stability and 
predicts no dependence on system size. Using the 
frictionless model and performing stability analy- 
sis within the pile, we reproduce the dependence 
of the stability angle on system size, particle size, 
and surface tension observed in our experiments. 
Additionally, we account for past discrepancies in 
experimental reports by showing that sidewalls 
can significantly increase the stability of granular 
material. 

The experimental apparatus consists of a clear plexi- 
glass drum which is rotated about an horizontal axis. We 
primarily use a drum with a diameter D of 28.5 cm and 
a width W which can be varied from Ocm to 14.5 cm. A 
drum with D = 12.5cm and W — 11.5cm is also used 
to vary system size. The rotation rate u> is varied from 
9.0 x 10 -4 rpm to 5.6 x 10 _1 rpm. Soda-lime glass spheres 
which have a density p of 2.4 gem' 3 , and with radius r = 
0.25 mm, 0.3 mm, 0.5 mm, and 1.5 mm and size dispcrsity 
within 0.1 mm were used. The drum was 40% filled with 
grains premixed with a small amount of liquid. We re- 
port the amount of liquid added in terms of the volume 
fraction, Vf, which is defined as the volume of the liq- 
uid divided by the volume occupied by the grains alone. 
The effect of surface tension of the liquid T was tested 
by using silicone oil and water which have r = 20 ± 1 
dyne cm -1 and T = 70 ± 1 dyne cm -1 , respectively. Sil- 
icone oil with viscosity v ranging from 5cS to 1000 cS 
is used to study its impact on the measured inclination 
angles. 

The drum is back lit so that light could pass through a 
thin layer of grains that tends to accumulate on the sides 
of the drum, but is unable to pass through the bulk of 
the pile. Images were acquired with a mega-pixel resolu- 
tion digital camera at a rate of three frames per second 
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FIG. 1: Avalanching regimes Glass spheres (radius r 0.5 
mm) mixed with silicone oil in a rotating drum (width W 
14.5 cm), (a) Stick-slip flow is observed for a rotation rate ui 
of 0.028 rpm and viscosity v = 5 cS. (b) Continuous flow is 
observed for rotation rate of 0.28 rpm. 



and used to determine the surface slope with automated 
code. At this rate, the slope of the pile changed no more 
than 0.2 degrees between frames. This error in measure- 
ment is not significant given that the slope of the pile at 
the moment of avalanche was distributed over about two 
degrees in a given run. At least 30 events were recorded 
to measure the mean values. 

A stick-slip avalanche regime or a continuous avalanche 
regime is observed depending on the rotation rate of the 
drum and the liquid viscosity [see Fig. IHa,b)]. In the 
stick-slip regime, the slope of the pile increases linearly 
with time until the heap reaches the maximum angle 9 m . 
At this point, the grains are observed to avalanche and 
the pile's angle decreases to the angle of repose, 6 r when 
avalanching stops. In the continuous avalanche regime, 
the grains were observed to flow continuously and the 
slope of the pile is approximately constant over time. By 
changing to and v, one can go from the stick-slip to con- 
tinuous avalanche regime. 9 m in the stick-slip regime is 
observed to be independent of lu and v within the fluc- 
tuations in the data and thus rate independent. On the 
other hand, 9 r and the surface angle in the continuous 
avalanche do depend on the viscosity, consistent with pre- 
vious observations . The overall behavior as a function 
to is similar to previous systematic investigations with a 
rotated drum [7J, but the magnitudes are significantly 
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FIG. 2: Impact of liquid volume fraction and drum 
width on stability (a) Measured surface angles as a function 
of Vf of silicone oil. (y = 5 cS, r = 0.5 mm, to = 0.028 rpm.) 
(b) 9 m is observed to decay exponentially as the width of the 
drum is increased. Error bar is indicated on one of the points. 
(y / = 8x 10~ 3 , a = 23.3 deg., b = 0.33cm- 1 , c = 30deg.) 



lower. 

We note that empirically, as function of Vf, 9 rn in- 
creases sharply, and saturates when Vf is on the order 
1 x 10 3 , 9 m [See figure E^a)]. The dependence of the an- 
gles on Vf is again qualitatively consistent with previous 
studies [a, 0, uut as noted earlier the overall increase 
and saturation is significantly lower. At very small Vf, 
9 m can be quite sensitive to Vf because of the inherent 
roughness of the grains and its impact on the shape of the 
liquid bridge 0. We performed our experiments in the 
saturation regime, so that such effects are eliminated, but 
we also chose a small enough Vf that the liquid does not 
drain to the bottom due to gravity, which would cause 
spatial inhomogeneity. 

We experimentally determined that higher angles can 
be attributed to significant wall effects when narrow 
widths are used as in most previous investigations. For 
example, the angles reported in Ref. are 15-20 degrees 
higher than those we report for similar r and T, but the 
apparatus was 3.2 cm wide. The observed dependence on 
W is plotted in Fig. H{b) an d is described by an exponen- 
tial fit. The effect of side walls on grains which are fully 
immersed in a liquid has been reported in Ref. @- The 
data in that case reaches the asymptotic value for a W of 
a few grain diameters because liquid bridges are absent. 
However, the side walls have an effect over many hun- 
dred times the diameter of the particle in the partially 
saturated case because of the added particle-particle cor- 
relations induced by the liquid bridges. In order to sim- 
plify the analysis, we examine the data only in the limit 
W > 11.5 cm where the side walls are unimportant. 




FIG. 3: Schematic of the liquid bridge model (a) A 

schematic of the system, note that L is defined as the length 
of the pile's surface, (b) a geometrical arrangement of a sphere 
resting on top of three spheres in contact used in the stabil- 
ity analysis. Liquid bridges located at the points of contact 
introduce forces along the segments joining the vertices abed 
of the tetrahedron at which the spheres are centered. 

8m of at least 90 degrees is predicted by Albert, et 
al. |2| if the ratio of the capillary to gravitational force, 
which is also called the bond number (Bo 0), is of order 
one or greater. Now the capillary force due to the liquid 
bridge bond is given by: 

F c = omTr (1) 

where, a is a dimensionless constant which depends on 
the size and shape of the liquid bridge between the parti- 
cles, and their separation 10illl| . Although there is some 
ambiguity on the value of a, it is of order 1, and therefore 
Bo is approximately 5 for glass beads with r = 0.5 mm 
mixed with silicone oil at Vf = 0.008. We do not observe 
9 m reach 90 degree; the piles typically avalanche when 
9 exceeds 30 degrees as can be noted from Fig. ^ We 
believe that the discrepancy arises from the assumption 
that a wet granular pile fails at the surface, as was also 
noted previously in Ref. H. 

Following Albert, et al. |2|, let us consider the stability 
of a sphere placed on top of three base spheres which are 
all in contact with each other (see Fig.[3J). In the absence 
of liquid bridges, the sphere becomes unstable when grav- 
itation force is outside the triangle formed by the base 
spheres. By averaging over all possible orientation an- 
gle 4> of the base triangle between and 60 degree, they 
showed that the average inclination 9d (the geometric, 
or dry angle of stability) is 23.8 degrees. At this angle, 
approximately half of the triangular bases at the surface 
can support a sphere and thus the pile will be stable. 
Measurements with dry glass, plastic or steel spheres in 
a rotated drum all show 9 m which are within a few de- 
grees of their prediction, and thus it appears that the 
complications associated with random packing and fric- 
tion between particles need not be taken into account to 
obtain reasonable results. 

Next, let us again consider a sphere resting on three 
base spheres as in Fig.^b) which is located at a distance 
x and angle 9 from the bottom of the inclined surface. For 
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FIG. 4: Comparison of data with liquid bridge model 

(a) The maximum angle of stability glass spheres mixed with 
water. [L — 27 cm). The solid line shows the fit to eq. (0 
with ao — 1.2. (b) The experimentally measured 9 — 8d are 
plotted against the theoretical values calculated from eq. 10. 
Key: o - silicone oil in the large drum, □ - water in the large 
drum, A - silicone oil in the small drum, V - water in the 
small drum. Filled shapes, open shapes, and shapes with a x 
correspond to r = 1.5 mm, 0.5 mm, and 0.3 mm, respectively. 



9 > 6d, the sphere will be gravitationally unstable. Now 
we assume that on an average, it is the component of the 
liquid bridge bond ab directed up the inclined plane tilted 
at 9 which is responsible for offsetting the unbalanced 
gravitation component down the plane. (We neglect the 
contributions of the bonds ac and ad because they are 
mostly directed towards the axis around which the top 
sphere rotates when it becomes unstable.) 

Because abed forms a tetrahedron, the projection of the 
force F c up the plane of the base triangle due to the liquid 
bridge between ab can be easily determined. The average 
component of the liquid bond force corresponding to the 
average orientation angle <fi of the base triangle is given 
by F c / ' \f2itaa.9d which is approximately 0.46F C . There- 
fore the shear stress that can be supported by the liquid 
bridges along the 9 plane is given by the force per sphere 
calculated times the average number of spheres per unit 
area. Now, the number of particles per unit volume is 
3/ p /47rr 3 , where f p is the packing fraction of the grains 
which is approximately 0.64 for spherical grains. There- 
fore the number of particles per unit area is equal to the 
number of particles per unit volume to the | power. We 
then have the shear stress as: 



( 3/p 



\/24tan0 d \4TTr 3 
where, F c is given by eq. QJ. 
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balanced weight of a thin vertical volume element above 
a plane tilted at 9 and located at a distance x from the 
bottom of the surface as shown in Fig.|3Ja). We can write 
the unbalanced weight to be given by: 



Mg sin {9~9 d ), 



(3) 



where, M is the mass of the volume element which has a 
length dx, width W, and a height xsm(6 m — 9)/ cos 8d- 
We take (9—9d) in eq. because the pile is geometrically 
stable up to 9d- Thus we find that the shear stress on 
the 9 plane from the weight of the volume element is 



xpgf p sin (9 m - 9) sin (9 - 9 d )/ cos Q d . 



(4) 



To simplify further discussion, we take the small angle 
approximation, and obtain the unbalanced stress to be 



xpgf P {9 m - 9)(9 - 9 d )/ cos9 d - 



(5) 



Now for each x it can be shown that this shear stress is 
maximum for 6 = m f 4 ■ This is simply because as one 
goes above 9d on one hand there is less materials above, 
but on the other hand the volume element is tilted further 
off equilibrium. 

Next, it is not clear if the pile will fail when the shear 
stress from the weight first exceeds the stress provided by 
the liquid bridges, or when the average shear stress due to 
the weight of the pile exceeds the stress provided by the 
liquid bridges. In reality, the correct physical situation is 
most likely combination of these two. Therefore, we set 
x — PL in eq. JSJ when balancing the stresses, where (3 
is a dimensionless constant between 0.5 and 1. j3 = 1, if 
adjacent volume elements cannot offset any stress by a 
particular volume element, and [3 = 0.5, if the stress is 
evenly distributed along the slip plane. Given the vari- 
ance seen in the stick-slip events [see Fig. ^a)], and the 
heterogeneity of pile at the granular level, [3 may vary 
from event to event. 

Combining these facts and balancing the two stress 
components, we obtain the condition for equilibrium as 



9 d 



( 9tt 



and, setting ao 



a/p, we obtain 



a cos \ 
24 tan 6 d pgr/3L J 



(6) 




ao cos 
V6 tan 9dpgrL 



(7) 



Now let us calculate the shear stress due to the un- 



Figure^Ja) shows a comparison of the measured 9 m as 
a function of particle size r and the fit to eq. J7J) after 
conversion of units to degrees. The data is well described 
by using ao = 1-2 ± 0.05. To further test the model, we 
measured 9 m in two different sized systems (L = 27 cm 
and L = 12.5 cm) with various combinations of Y and r. 
Fig- Elk) shows a plot of the experimentally measured 9 m 
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above 9d compared with eq. J7J), with ccq — 1.2. Excellent 
overall agreement is observed. 

We note that the small angle approximation used to 
simplify eq. is reasonable over the angles tested. The 
analysis can be extended to higher angles without us- 
ing this approximation, but 9 m has to be obtained by a 
numerical solution of the corresponding equation. In de- 
veloping our analysis, we have ignored the contribution 
of friction which will also help stabilize the wet pile, but 
just as in the dry case, we may expect the friction con- 
tribution to be small as a grain rolls as it gets dislodged 
from the triangular base. A more detailed knowledge of 
how the failure develops may clarify this issue. 

In conclusion, our analysis captures the observed sta- 
bility dependence on grain size, system size, and surface 
tension. It would be of great interest to consider how well 
this analysis extends to piles composed of non-spherical 
and multi-sized grains as in natural sand. 
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